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In order to fulfill the recent thrust of Astatine radionuclides in the field of nuclear medicine 
various production routes have been explored in the present work. The possible production routes 
of 2°»-2"At comprise both light and heavy ion induced reactions at the bombarding energy range 
starting from threshold to maximum 100 MeV energy. For this purpose, we have used the nuclear 
reaction model codes TALYS, ALICE91 and PACE-II. Excitation functions of those radionuclides, 
produced through various production routes, have been calculated using nuclear reaction model 
codes and compared with the available measured data. Contribution of various reaction mechanisms, 
like, direct, preequilibrium and equilibrium reactions, to the total reaction cross section has been 
studied using the codes. Result shows that equilibrium reaction mechanism dominates in all cases 
over other reaction mechanisms. 

PACS numbers: 24.10.-i, 25.55.-c, 25.70.-z 



I. INTRODUCTION 

In the last decade semi-metallic radioactive element astatine (At) has drawn ample attention 
in the field of nuclear medicine. Nuclear physics has an important role in developing this field 
as nuclear reaction is the method of choice to produce At radionuclides. Most of the identified 
radioisotopes of At are short lived. ^-'^^At, which has half life, Ti/2=8.1 h, is the most stable isotope 
among others. Though At radioisotopes are produced artificially, however, in the Earth's crust, 
small amount of At may last in equilibrium with uranium and thorium. At was discovered about 
80 years ago, but, studies on its chemical properties faced practical constraint due to its scarcity 
and short half life. Some recent reports on chamemical separation of At radionuclides from the 
bulk, which is also an important step followed before application, has been cited in ref 
list of At radionuclides, which may find possible use in nuclear medicine, having suitable half life 
is presented in Table [J 

Among others, ^^^At is most promising radionuchde for targeted therapy due to its suitable half 
life and relatively high linear energy transfer (LET) of a-particles in human tissues. ^^^At has two 
major decay modes both of which eventually emits high energy a-particles suitable for targetted 



TABLE I: The nuclear properties [4| of At radionuclides useful in the field of nuclear medicine 
Radionuclide T1/2 Decay mode Q(MeV) 7(keV)(%) 
^"At 7.214 h e 58.2%, a 41.8% 5.867 687(0.246) 

At 8.1 h e 99.82%, a 0.18% 5.524,5.361 1181.43(99.3),245.35(79), 1483.35(46.5) 
At 5.41 h e 95.9%, a 4.1% 5.647 545.03(91),781.89(83.5), 790.2(63.5) 
20«At 1.63 h e 99.45%, a 0.55% 5.641 685.2(97.89),660.1(90), 177.0(46) 
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therapy. The a-decay mode emits a-particle of energy 5.867 MeV leaving ^°^Bi (Ti/2=32.9 y) as 
residual, which decays to stable ^"^Pb. Other decay mode produces ultra short lived daughter, 
2ii9po (T^^2 = 516 ms, J =9/2+), which decays to stable ^o^pb emitting 7.45 MeV a-particle along 
with the emission of X-rays. The metastable state of 2ii"ipo (T1/2 — 25.2 s, J =25/2+) is not 
populated by ^^^At because of its high spin and low transition probabilities. The 7-emitting ra- 
dionuclides 209, 210^^ have mostly e-capture decay mode and produce long lived a-emitting 209,2iopQ 
daughter. These At radionuclides may find uses in diagnosys due to their suitable half life and 
high intensity 7-rays. 

Radionuclidic purity is an important issue in nuclear medicine. Therefore, At radionuclide for 
medical application has to be obtained in its purest form to avoid unnecessary radiation exposure 
from impurities. In general, production route and the chemical separation technique are the two 
major sources of impurity (radionuclides as well as stable) that is propagated to the particular 
isotope of interest during its production and purification procedures. Production route hosts ra- 
dionuclidic impurity and their corresponding decay products including the decay of compound 
nucleus. Optimization of nuclear reaction parameters, like, suitable target-prijectile combination, 
bombarding energy of projectile, thickness of target, etc., is a probable way to reduce certain pro- 
duction of impurities. For exarnple, the most commonly used production route of At radionuclides 
is 209Bi(Q;,xn)209-2iiAt 0, 0, 0, B S, 0, [nl. These At isotopes and their daughter products 



(radioactive and stable) can be considered as impurity when aimed for a particular isotope. Long 
lived impurities are hazardous to human health if the amount goes above a certain level while 
the stable heavy metals, Pb, Tl, which are the end products of At radioisotopes, are sometimes 
carcinogenic and affects human metabolism. 

At radionuclides can be produced directly by nuclear reactions using suitable target-projectile 
combination or from the decay of its precursor Rn isotopes. Light and heavy ion induced reactions 
may be used for the production of At radionuclides in low and medium energy accelerators. High 
energy proton induced spallation reactions on heavy targets like, ^^^U, ^'^^Th can also be used to 
produce ^o^-^HAt. In the indirect method, ^o^-^iiRn (Ti/2=28.5 m, 2.4 h and 14.6 h respectively) 
isotopes are produced either by light charged particle induced reactions or by spallation reaction 
induced by high energy proton. At isotopes are produced as a daughter product of Rn isotopes. 
Various production routes of 209-211^^ escort different type of impurities whose quantification is 
important for their application in nuclear medicine. 

The present work is to explore the possible production routes of 209-211^^ radionuclides from 
both light and heavy ion induced reactions. Production of At radionuclides is preferred from a- 
induced reaction on ^^^Bi because of high yield required for theraputic purpose. However, produc- 
tions via light heavy ion induced reactions may serve some other scientific purposes, like, biological 
studies which requires minute amount of radioactivity. We aim in this paper to investigate theo- 
retically the various possible production routes of ^^^"^^^ At from both light and heavy ion induced 
reactions and the reaction mechanisms involved in each route. A comparison of theoretically pre- 
dicted production cross sections has been made with the reported experimental values wherever 
available. 

Production of 209-211^^ radionuclides from light {a and '^He) and heavy ion (^^^Li, ^Be) induced 
reactions on heavy natural targets like Bi, Pb, Tl, involves mostly preequilibrium (PEQ) and 
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equilibrium (EQ) reaction mechanisms in the low and medium energy range. In case of high 
projectile energy, direct reaction (DIR) mechanism may also have some contribution. A composite 
nucleus is formed when a target is bombarded by the projectile of certain energy. Particle emission 
may occur from any stage of the relaxation process of the target-projectile composite system leaving 
the residual in an excited state. The excited residual nucleus is cooled either by particle emission, if 
the residual excitation is higher than the particle emission threshold giving birth to a new residual 
nucleus, or by gamma ray emission. 

In case of direct reaction, particles are emitted due to large momentum transfer from the very 
first interaction between the projectile and the target. As a result ejectile carries high energy 
leaving the residual at lower excitation. When incoming energy of the projectile is shared by 
all the nucleons of the composite system, equilibrium condition is reached. Particles are emitted 
from the compound nucleus due to statistical fluctuation in energy. Therefore, compound nuclear 
emissions are low energy emissions and the residual is left at high excitation energy. 

Between the two extremes, there is an intermediate stage of reaction which exhibits both direct 
and compound like features. These reactions are referred to as PEQ processes. PEQ emissions 
take place after the first stage of the reaction (DIR) but long before statistical equilibrium of 
the compound nucleus is attained. PEQ processes cover a significant part of the reaction cross 
section for incident energies between 10 and 200 MeV. The production cross section of the desired 
radionuclide is the resultant of all the reaction mechanisms involved. However, contribution of 
direct reaction mechanism is comparatively small in the energy range up to 100 MeV considered 
for the production of At-radionuclides in the present work. 

Light and heavy ion induced nuclear reactions have been studied using the nuclear reaction model 
codes TALYS Q, ALICE91 0, Q and PACE-II [li. The code TALYS takes into account 
all three types of reaction mechanisms; direct (Coupled channel analysis. Distorted Wave Born 
Approximation (DWBA) etc.), preequilibrium (one and two component exciton model) and equi- 
librium (Hauser-Feshbach model) , PACE-II accounts only evaporation reactions based on Hauser- 
Feshbach model [l3| and ALICEQlconsiders preequihbrium(hybrid model) along with the equi- 
librium reactions (Weisskopf-Ewing model) [l7[. A list of reactions studied in the present work 
is presented in Table [Til Use of three codes enables us to understand the reliability as well as 
predictablity of the codes. 

In the next section we give a brief description of the codes TALYS, ALICE91 and PACE-II and 
the input option used for the present calculation. Sec III shows the remarkable features of the 
codes. In sec IV we discuss the results obtained in the present study. 



II. BRIEF DESCRIPTION OF NUCLEAR REACTION MODEL CODES 

A. ALICE91 

We have used the code ALICE91 0, [l3| to calculate the excitation function of At radionuclides 
produced through various reaction mechanisms. Geometry dependent hybrid model [14l.ll8l.ll9l.l2C| 
has been used to calculate PEQ emissions and Weisskopf-Ewing formalism [l3| for EQ emissions. 
ALICE91 does not take care of DIR emissions. The hybrid model is the combination of exciton 
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TABLE II: List of nuclear reactions studied along with corresponding Q values and threshold energies 



Reaction Q value (MeV) Eth(MeV) 



2°''Bi(a, 2n)2"At 


-20.33 


20.72 


2°9Bi(a, 3n)2iOAt 


-28.07 


28.61 


2°''Bi(a, 4n)209At 


-35.24 


35.92 


2°''Bi(='He, n)2"At 


0.25 





2°^Bi(3He, 2n)2"'At 


-7.50 


7.61 


2°^Bi(3He, 3n)2iOAt 


-14.66 


14.87 


^o^PbC'Li, 3n)'"At 


-20.23 


20.81 


^"^PbC^Li, 4n)^^°At 


-27.98 


28.78 


^ospbC^Li, 5n)2"«At 


-35.14 


36.15 


208pb(^Li, 4n)2"At 


-27.48 


28.41 


208pb(^Li, 5n)2iOAt 


-35.23 


36.41 


208pb(^Li, 6n)209At 


-42.39 


43.82 


203Tl(9Be, n)2iiAt 


-10.84 


11.32 


2°3Tl(«Be, 2n)2i''At 


-18.58 


19.41 


203Tl(SBe, 3n)2iOAt 


-25.75 


26.89 


205Tl(SBe, 3n)2"At 


-25.04 


26.14 


205Tl(SBe, 4n)2i»At 


-32.79 


34.23 


205Tl(«Be, 5n)2»9At 


-39.95 


41.71 
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assumes that the target-projectile 



model 2lj and Boltzmann master equation approach 22 
composite system proceeds through two body interaction process. Each stage of the relaxation 
process is designated by the total number (n) of excited particles, sum of excited particles (p) and 
hole (h). In each two body interaction, p-h pair may be created or annihilated or redistribution 
of energy takes place without changing the number. Hybrid model uses "never come back" ap- 
proximation, i.e, the model assumes only p-h pair is created in each interaction. Hybrid model 
explicitly determines the pre-emission energy distribution of the excited particles which helps to 
estimate high energy emissions more accurately. Geometry dependent hybrid model includes the 
nuclear surface effects 



IJ, [191, l20[ . The PEQ emission cross section for a particular ejectile x with 



energy is given by 



2 oo n 

Here, A is the de-Broglie wave length of the projectile, Ti is the transmission coefficient of the 
Z*'' partial wave, Z?„ is the depletion factor of the n*'' exciton state, that is, probability of reaching 
n exciton state without prior emission and is number of x type excited nucleon present in 
it. The number no and n are the initial and equilibrium exciton numbers respectively. The ratio 
Nn{l, U, ex) /Nn{l, Ec) is the probability of finding x type nucleon in the n exciton state with energy 
(ex+Bx) where Bx is the separation energy of x. The factor ^ (^^"^^^^^^ ) is the emission probability 
of x with energy e^,. Xti^x) is the two-body interaction rate. The emission rate Xd^x) is calculated 



n N„(l,E^) 



Ac(e=,)+At(ex) 



(1) 



by 



H 



KK^x) = ^-3 i^j 



{2Sx + l)^.xex (ex) 

where, Sx is the intrinsic spin of x, fi^ is the reduced mass, ainv is the inverse cross section of the 
ejectile x with energy ex being absorbed by the residual and g is the single particle level density 
of the composite nucleus. Equilibrium emission cross section is calculated using Weisskopf-Ewing 
formalism as 

(3) 

cTcomp is the compound nuclear formation cross section, a is the level density parameter and U is 
the available excitation energy of the compound nucleus after the preequilibrium emissions, a comp 
is calculated as a comp — o'ats ~ ctpeq where aabs is the absorption cross section of the projectile 
in the target and apEQ is the total preequilibrium emission cross section. 



1. Input options used 
The calculation has been performed using the code ALICE91 0,0 

with geometry dependent 

hybrid model for PEQ emissions and Weisskopf-Ewing formalism for EQ emissions. Neutron, 
proton, alpha and deuteron emissions are considered from the residual nuclides of 12 mass unit 
wide and 10 charge unit deep including the composite nucleus. Fermi gas level density has been 
used for the calculation of reaction cross section. Reverse channel reaction cross sections have been 
calculated using the optical model. Default value of the level density parameter, A/9 has been 
chosen for the calculation. Total number of nucleons in the projectile has been chosen as the initial 
exciton number for the cross section calculation. 



B. TALYS 



The code TALYS 



12| has been used to calculate the production of At radionuclides through a- 
particle and ■^He-induced reactions on ^"^Bi target. We have adopted two-component exciton model 
for the estimation of PEQ emissions, detail Hauser-Feshbach formalism [igJ for the EQ emissions 
and Coupled Channel analysis for DIR reaction process. In case of PEQ reactions of projectiles and 
ejectiles heavier than nucleon (d, t, ^He and a-particles) direct like reaction processes; stripping, 
pickup and knockout reactions have important contribution in addition to the exciton model and 
are treated in the frame work of Kalbach's phenomenological model ^] as these processes are not 
included in the exciton model. 

In the two-component exciton model, total exciton number is divided into proton (pp) /neutron 
(pn) particle number and proton (hp) /neutron (ft,„) hole number. np{— pp + hp) and n„(= Pn + hn) 
are the proton and neutron exciton numbers respectively. PEQ cross section of the particle x with 
energy is calculated as 

Pp p„ 

C^PEQi^x) = Cr'~^^ D{Pp,hp,Pn,hn)W{pp,hp,Pn,hn,ex)T{pp,hp,Pn,K) (4) 

Pp=Pp Pn=Pn 
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where, a'~^^ is the composite nucleus formation cross section, D{pp, hp,pn, hn) is a fraction of 
PEQ population that has survived previous emission, W{pp, hp,pn, hn, ex) is the emission rate of x 
with energy and t is the mean life time of that state. Compound nuclear emission is calculated 
in TALYS by two different mechanisms; binary compound cross section that is the capture of 
the projectile in the target nucleus subsequently emitting particle or 7-ray at low energy and 
multiple compound emissions, that is, multiple emissions of highly excited residual nuclei formed 
after the binary reaction. Width fluctuation is taken into account in case of binary compound 
cross section. Multiple compound and precompound emissions are considered in the framework 
of Harser-Feshbach and exciton model. TALYS contains both macroscopic and microscopic level 
densities. Macroscopic level densities are 

(i) Constant Temperature model [2^ and Fermi gas model 

(ii) Back-shifted Fermi gas model and 

(iii) generalized superfluid model. 



1. Input options used 

The code TALYS uses detail Hauser-Feshbach formalism for multiple particle emission in the 



evaporation calculation along with width fluctuation corrections using Moldauer model [2^ [27[ . 
The option for multiple pre-equilibrium emissions has been chosen within the frame work of two- 
component exciton model calculation with surface correction to account the pre-equilibrium emis- 
sions. We have calculated level densities using a combination of the Constant Temperature Model 
by Gilbert and Cameron [2^ and the Fermi gas model. The total excitation energy range is divided 
into two regions: low energy part and high energy part. Low energy part starts from MeV to the 
certain energy, Emid, up to which constant temperature law is valid. The high energy part starts 
above the Emid and Fermi gas model is used to calculate level densities in this zone. 



C. PACE-II 

The code PACE-II is the modified version of the Monte-Carlo code Projection Angular- 
momentum Coupling Evaporation. The deexcitation process of the excited nuclei is calculated 



using the modified version of the code JULIAN [28|, which follows the correct procedure of an- 
gular momentum coupling at each stage of deexcitation. The code can be used to study the high 
energy heavy ion reactions. The transmission coefficients for light particle emission are determined 
from the optical model potential where all the optical model parameters are kept as default. The 
shift in coulomb barrier during deexcitation is accounted by calculating the transmission coefh- 
cients at an effective energy determined by the shift. The code internally decides level densities 
and masses it needs during deexcitation. Gilber-Cameron level density is taken as default. Fission 
is considered as a decay mode where the fission barrier can be changed accordingly in the program. 
The default fission barrier is the Cohen-Plasil-Swiatecki rotating liquid drop fission barrier. Saddle 
point level density is proportional to the a parameter. For compound nucleus reactions distri- 
butions are determined as function of angle around recoil axis and for deep inelastic products. 
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angular distributions are determined as functions of the angle axis perpendicular to the reaction 
plane. Compound nuclear fusion cross section is determined by using the Bass method 29]. A non- 
statistical yrast cascade 7-decay chain has been artificially incorporated to simulate 7 multiplicity 
and energy results. 



III. IMPORTANT FEATURES OF THE CODES 

(1) TALYS uses exciton model to calculate PEQ emissions followed by Hauser-Feshbach formal- 
ism for EQ emissions. ALICE91 uses hybrid model for PEQ and Weisskopf-Ewing formalism for 
EQ emissions 



31 



(2) TALYS relies on DIR emissions treated by coupled channel analysis, DWBA 30|, Giant reso 
nance 



32l |. etc. However, ALICE91 does not consider DIR emissions. 

(3) Direct like processes, for example, stripping, pickup, knockout, etc., have important role in the 
nuclear reactions involving cluster (d, t, ^Ke and a-particles) as projectile and ejectile. TALYS 
includes these processes in the frame work of Kalbach's phenomenological model. ALICE91 does 
not take into account the above described processes. 

(4) TALYS considers n, p, d, t, '^He and a-particle as both projectile and ejectile with single and 
simultaneous emission of any order in PEQ and EQ emissions. ALICE91 considers light particles 
up to a as projectile but only n, p, d and a-particle are considered as ejectiles. Single and si- 
multaneous emission of two nucleons are taken into account in PEQ emissions and in case of EQ 
reactions, n, p, d and a-particle emissions are calculated from successive residuals. PEQ emissions 
from heavy ion induced reactions can also be calculated using hybrid model in ALICE91with nec- 
essary modifications. 

(5) PACE-II is a Monte-Carlo statistical code based on Hauser-Feshbach formalism and the decay 
channels consider statistical emission of 7-rays, n, p and a-particles. 



IV. RESULTS AND DISCUSSIONS 

In the present report excitation function of 209-2ii^|- j^g^yg been calculated using the reaction 
model codes TALYS, ALICE91 and PACE-II produced through various possible target projectile 
combination in the incident energy range starting from threshold to 100 MeV. Table HI represents 
the details of reactions with the reaction threshold values. Production of ^°^~^^^At through ^°^Bi(a, 
xn) and ^°^Bi(^He, xn) reactions have been studied using the code TALYS and ALICE91 and 
have compared with the available experimental cross section. Other possible production roots of 
At radionuclides via heavy ion induced reactions have been studied using the codes ALICE91 and 
PACE-II and have also compared with the available measured data. It has been seen that all the 
measured excitation functions of -At flHHH produced through -^Bi(a, 2n)-^At reaction 



are in good agreement with each other. 

It has been seen earlier that a-particle may interact with the target nucleons in different ways. 
Among others, 95% contribution comes from the two modes of interaction to the total reaction cross 
section; (i) a-particle dissolves into the four constituent nucleons in the nuclear force field (no=4; 
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Incident energy, E [MeV] 



FIG. 1: Comparison between the calculated excitation functions of ^^'^ At-radionuclides from the reac- 
tion ^°'^Bi(a, xn)^"^"^"At using the codes TALYS (lines) and ALICE91(lines with symbols). 



4p-0/i) Sy, |37| (ii) a-particle as a whole interacts with a target nucleon and creates a particle- 



hole pair that is the starting point of the PEQ cascade (no=6; 5p-lh). It has also been found 
that both these initial configurations have proportional contribution to the angular distribution 
of the ejectiles 38|. ALICE91 has option to select initial exciton number, no in the input. We 
have performed separate calculation of excitation function for both the options hq—A and 6. The 
resultant excitation functions of ^°^~^^^At have been calculated taking the 70% contribution from 
no=6 and 30% of that from no=4 according to the observation of Gadioli et al 0. 

Figure [T] shows a comparison between the formation cross section of ^°^~^^^At produced from 
209Bi via (a, 2n), (a, in) and (a, 4n) reaction channels calculated using TALYS and ALICE91. In 
the TALYS calculation, optical model has been used to calculate inverse reaction cross section with 
default level density option, which is the combination of Constant Temperature model of Gilbert 
and Cameron [25] and Fermi gas model. In the ALICE91 calculation we have used optical model 
with Fermi gas level density option. All three isotopes of At are well separated from each other and 
the maximum production of At and ^^°At have also been seen at the at the incident energies, 30 
MeV and 40 MeV respectively in both the results. But, ALICE91 produces 20% more ^^^At and 
10% less At than TALYS. Maximum production of ^"^At has been seen at 45 MeV in ALICE91 
and that at 50 MeV in TALYS result. 

To investigate the results obtained from the two model codes, we have studied the contributions 
from the reaction mechanisms involved in a-induced reaction on ^^^Bi target. Contributions from 
PEQ, EQ and DIR reactions to the total reaction cross section have been shown against incident 
projectile energy from TALYS and ALICE separately in Figure [2] and [H In figure [21 we see 
that there is almost no contribution from direct reaction for incident energy up to 100 MeV. 
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Incident energy, E [MeV] 



FIG. 2: Contributions of different reaction mechanisms to the total reaction cross section as obtained from 
TALYS. 




Incident energy, E [MeV] 

FIG. 3: Contribution of total PEQ and EQ reactions to the total cross section as obtained from ALICE91 
for two initial exciton numbers (no=4; 4p-0h and no=6; 5p-lh) in the a-induced reaction. 
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No PEQ reaction contribution has been seen in low incident energies, however, PEQ component 
increases with the increasing incident energy. The most dominating component is the EQ reaction 
in the range considered for the present work. It foUows similar increasing trend up to 50 MeV 
with the total reaction cross section, which shows increasing trend through the entire energy range 
considered. The rate of increment of EQ cross section is slow at the higher incident energies as PEQ 
component is increasing with the increasing incident energy. In case of ALICE91 in figure [31 we see 
that EQ contribution is significant only at low energies for both the initial configurations {Ap-Oh 
& 5p-lh) of a- induced reaction. PEQ component is increasing rapidly with the increasing incident 
energy for both the initial exciton numbers. The total cross section also shows the increasing 
trend throughout the range. PEQ contribution from 5p-lh configuration is less than that of Ap-Oh 
configuration of a up to the energy range 60 MeV. This is due to the fact that in case of 5p-lh 
configuration, PEQ emission starts in the hybrid model calculation from no=6, which is one stage 
later that no=4. As a result, PEQ contribution is larger for no=4 than the other. From the plot 
we see that total PEQ cross section does not increase above 60 MeV for both the cases. This is 
because, in the code ALICE91 has a cut off value {— no+18) if it exceeds this value. The case of 
EQ reaction is just reverse between no=4 and no=6. This is because after sufficient preequilibrium 
emission from the Ap-Oh configuration total available excitation becomes low which in turn reduces 
the EQ emission. 

Though the favored production route of ^"^"^^"'^At is from a-induced reaction on Bi target, it 
is also possible to produce 209-211^^ using light heavy ion induced reaction. In this report, we 
have found possible production routes of 209-211^^ from ^'^Li and ^Be induced reactions on "'^*Pb 
and "'^'Tl targets. Figure [4][6] represent the comparative study of the formation cross section of 
209-211^^ produced through (''Li, xn), (^Li, xn) and (^Be, xn) reactions on "'»*Pb and "'»*T1 targets 
using the codes ALICE91 and PACE-II. 

The four naturally occurring isotopes of Pb are ^ospb (52.4%), 207pb (22.1%), ^°^Ph (24.1%), 
204pb (1.4%). Tl has two naturally occuring isotopes, ^ostI (70.5%) and ^03x1 (29.5%). The 
formation cross sections have been calculated separately from each isotope of Pb using ^Li and ^Li 
projectiles and that from ^Be induced reactions on Tl isotopes. The total formation cross section is 
calculated taking the weighted average of all the naturally occulting isotopes. The same procedure 
has been followed for the results obtained from both the codes ALICE91 and PACE-II. 

The total formation cross section obtained from ALICE91 is the sum of PEQ and EQ cross 
sections. In case of PEQ calculation in the framework of hybrid model for heavy ion induced 
reaction we have considered the initial exciton number as the total number of the nucleons in 
the projectile. The EQ cross section is calculated using the Weisskopf-Ewing formalism. The 
formation cross sections obtained from PACE-II is only from the evaporation calculation following 
the Hauser-Feshbach formalism. 

In figure [4] and [5l production cross sections of 209-211^^ obtained from ALICE91 are larger than 
that obtained from PACE-II but the maximum production cross section is obtained at a particular 
incident energy for both the model calculations. Production of ^^^At is not significant in the ^Li- 
induced reaction on ""'Pb target. However, 209-2io^-(. ^j^l be produced in good amount in both 
''Li and ^Li induced reaction on natural Pb target. Figure [6] shows the production cross section 
of 209-211^^ radionuclides from ^Be induced reaction on '"^*T1 target. Almost no production of 
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FIG. 4: Comparison between the formation cross sections of ' 
target using ALICE91 and PACE-II. 



"^At from ^Li induced reaction on ""^Pb 



^^^At is seen in this case with sufficiently good production of two other radionuchdes. In all these 
three cases (figure |3][ni), ALICE91 prediction is larger than PACE-II. This may be due to the added 
contribution from PEQ reactions in ALICE91. 

In fi gur e [7]|51 we have compared the measured production cross sections of 2" At lO] and ^lOAt 



10 



4C 



4l[ produced from a-induced reactions on ^^^Bi target through (a, 2n) and {a, 3n) channels 



reported by several authors with the theoretical predictions obtained from TALYS, ALICE91 and 
PACE-II. Figure n.h„w, .ha. TALYS pertec.i; matches with the measured data of -A. 

in the range shown. ALICE91 also matches with the measured data except in the low energies. 
PACE-II reproduces the data up to 32 MeV. It underpredicts the data showing a abrupt fall in 
the cross section above 32 MeV. This comparison reveals that PEQ reaction mechanism plays an 
important role in reproducing the measured data. We have seen in figure [2] that in the energy range 
up to 100 MeV, DIR reaction has absolutely no role. As TALYS and ALICE91 account the PEQ 
reactions along with EQ reaction, they produce the measured data satisfactorily. Figure [8] shows 
that TALYS also evaluates the measured formation cross section of ^^^At. ALICE91 produces the 
data well but not below 32 MeV. PACE-II shows the same nature as the production of ^^^At. 

Figure [9] compares the measured excitation function of -^^^At from •^He- induced reaction 42ll43, 
44| on the ^^^Bi target with the TALYS predicted value. The trend of the measured data is well 



reproduced by TALYS, but its absolute value underpredicts the measured data by the factor of 
five. ALICE91 absolutely fails to produce the data. Figure [TUlfTT] present comparison between 
measured data of ^^At and ^lOAt produced from ^°^Ph{^U,3nf^^ At and 208pb(6Li,4n)2iOAt 



45j reactions with the theoretical values from ALICE91 and PACE-II. From figure [TOl we observe 
that formation cross section of ■^^^ At is well reproduced by the calculated values of ALICE91 within 
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FIG. 7: Comparison of measured formation cross section jiot] of ^^^At from a-induced reaction on ^"^Bi 
with the TALYS, ALICE91 and PACE-II predictions. 

the uncertainty range. But PACE-II only reproduces the higher energy side of the spectrum. In 
case of formation cross section of ^^°At in figure [TTl ALICE91 and PACE-II both reproduce the 
trend of the measured data. Absolute values of cross section obtained from PACE-II produce the 
measured data, whereas ALICE91 slightly overpredicts the data. 



In a similar fashion, we have compared production of ^^^At [46| and ^^°At [46] from ^Li-induced 
reaction on •^"^Pb target. The trend of measured cross section of ^^-'^At in figure [12] is well repro- 
duced by ALICE91 but absolute value of cross sections over predicts the measured data below 45 
MeV. PACE-II neither produces trend or the absolute values of cross section except for a particular 
energy 45 MeV. In case of figure [TH we see that ALICE91 explains the measured data well but 
PACE-II under predicts the data mostly. 



The study of the heavy ion induced production of 209-211^^ radionuclides, we can conclude that 
PEQ reaction mechanism plays an important role in case of light heavy ion induced reactions in 
the energy range studied and performance of the ALICE91 code is better in most cases. PACE- 
II mostly underpredicts the measured at high energy side where PEQ reaction contributes. As 
the code does not consider the PEQ reaction in its calculation, underpredicts the measured data. 
TALYS and ALICE91 calculated excitation functions of At radionuclides from a-induced reaction 
are comparable and both agree with the experimental data at higher energies, but the cross sections 
at low energies are well reproduced by TALYS. Heavy ion induced production cross sections of 
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209-211^^ obtained from ALICE91 are mostly in good agreement with the measured data, whereas 
PACE-II calculations underpredict them. This may be due to the significant contribution from 
PEQ reactions which is accounted in ALICE91. The study on the contribution of various reaction 
mechanisms to the total cross section agrees with fact. Direct reaction has no role in the production 
of At radionuclides from both light as well as heavy ion reactions. 



This work has been carried out as a part of the Saha Institute of Nuclear Physics-Departmet 
of Atomic Energy, XI five year plan project "Trace Analysis: Detection, Dynamics and Spcciation 
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